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For more than three decades, the inner crust of neutron stars, formed of a solid lattice 
Q , of nuclear clusters coexisting with a gas of electrons and neutrons, has been tradition- 

l_l 1 ally studied in the Wigner-Seitz approximation. The validity of this approximation is 

^ ' discussed in the general framework of the band theory of solids, which has been recently 

j^j i applied to the nuclear context. Using this novel approach, it is shown that the unbound 

neutrons move in the crust as if their mass was increased. 
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1. Introduction 

cn _ 

, The interpretation of many observational neutron star phenomena is intimately 

related to the properties of neutron star crusts: 1 pulsar glitches, giant flares and 
toroidal oscillations in Soft Gamma Repeaters, X-ray bursts and superbursts, ini- 
tial cooling in quasi-persistent Soft X-ray Transients, precession, gravitational wave 
emission. Besides the ejection of neutron star crust matter into the interstellar 
medium has been invoked as a promising site for r-process nucleosynthesis, which 
still remains one of the major mysteries of nuclear astrophysics. 2 

The outer layers of the crust, at densities above ~ 10 6 g.cm -3 , are formed of 
a solid Coulomb lattice of fully ionized atomic nuclei. Below ~ 6 x 10 10 g.cm -3 , 
the structure of the crust is well established and is completely determined by the 
experimental atomic masses. 3 It is found that nuclei become increasingly neutron 
rich with increasing depth. At densities above ~ 6 x 10 10 g.cm -3 , nuclei are so 
exotic that experimental data are lacking. Nevertheless, the masses of these nuclei 
could be determined in the near future by improving experimental techniques. At 
density ~ 4.10 11 g.cm -3 , calculations predict that neutrons "drip" out of nuclei 
forming a neutron liquid. Such an environment is unique and cannot be studied in 
the laboratory. We will refer to the nuclei in the inner crust as nuclear "clusters" 
since they do not exist in vacuum. The inner crust extends from ~ 4.10 11 g.cm -3 
up to ~ 10 g.cm -3 , at which point the clusters dissolve into a uniform mixture 
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Fig. 1. In the Wigner-Scitz approximation, the crust (represented here as a two dimensional 
hexagonal lattice) is divided into independent identical spheres, centered around each site of the 
lattice. The radius of the sphere is chosen so that the volume of the sphere is equal to where 
71n is the density of lattice sites (clusters). 

of electrons and nucleons. Near the crust-core interface, the clusters might adopt 
unusual shapes such as slabs, rods or more complicated structures. 1 ' 4 

The structure of the inner crust has been studied using different approximations 
and nuclear models. The current state-of-the-art is represented by self-consistent 
quantum calculations which were pioneered by Negele&Vautherin in 1973. 5 In these 
calculations, the crust is decomposed into an arrangement of identical spheres cen- 
tered around each cluster as illustrated in Fig. [1] Each sphere can thus be seen as 
one big "nucleus" so that the usual techniques of nuclear physics can be directly 
applied to study neutron star crust. However, it has been recently found that this 
approximation, proposed a long time ago by Wigner-Seitz 6 in solid state physics, 
leads to unreliable predictions of the structure and composition of the crust es- 
pecially in the bottom layers. 7 Besides this approach does not allow the study of 
transport properties which are essential for interpreting observations. A more elab- 
orate treatment of neutron star crust beyond the Wigner-Seitz approximation is 
therefore required. 8 

2. Band theory of solids 

It is usually assumed that neutron stars are composed of cold catalyzed matter 
(for a discussion of deviations from this idealized situation, see for instance Ref. 1). 
Besides, it is rather well-established that the ground state of dense matter (at zero 
temperature) below saturation density possesses the symmetry of a perfect crystal. 

A crystal lattice can be partitioned into identical primitive cells, each of which 
contains exactly one lattice site. The specification of the primitive cell is not unique. 
A particularly useful choice is the Wigner-Seitz or Voronoi cell defined by the set 
of points that are closer to a given lattice site than to any other. This cell reflects 
the local symmetry of the crystal. It follows from the translational symmetry that 
a primitive cell contains all the information about the system. The single particle 
states of any species q in a periodic system are characterized by a wave vector A;. 
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If the wavefunction tpi? (r) is known inside one cell, the wavefunction in any other 
cell can be deduced from the Floquet-Bloch theorem 9 

^(r + T)=e ik - T ^\r), (1) 

where T is the corresponding lattice translation vector (which translates the initial 
cell to the other). Apart from the wave vector A;, the single particle states are labelled 
by a discrete index a (principal quantum number) so that the single particle energy 
spectrum consists in a series of "bands" or sheets in fe-space. This band index 
accounts for the local rotational symmetry of the lattice. 

It can be shown that the single particle states (therefore the single particle 
energies) are periodic in A;-space 

viSjc(r) = ^' ) (r), (2) 

where the reciprocal vectors K are defined by 

K T = 2ttN , (3) 

./V being any positive or negative integer. The set of all possible reciprocal vectors 
define a reciprocal lattice in fc-space. Equation @ entails that only the wave vectors 
k lying inside the so-called first Brillouin zone (i.e. Wigner-Seitz cell of the reciprocal 
lattice) are relevant. The first Brillouin zone can be further divided into irreducible 
domains by considering the rotational symmetry of the lattice. 

The inner crust of neutron stars is constituted of neutrons, protons and elec- 
trons. Unlike the situation in ordinary solids (ordinary meaning under terrestrial 
conditions), the electronic properties in neutron star crust are very simple. The 
matter density is so high that the Coulomb energy of the electrons is negligible 
compared to their kinetic energy. The electrons can thus be treated as a degenerate 
relativistic Fermi gas. 4 Since all the protons are bound inside nuclear clusters (ex- 
cept possibly in the bottom layers of the crust), the proton single particle states are 
almost independent of k and to a good approximation can thus be described only 
by the discrete principal quantum number a like in isolated nuclei. In contrast, the 
effects of the nuclear lattice on the neutrons (i.e. dependence of the states on both 
k and a) have to be taken into account because some neutrons are unbound. In 
the following, we will consider the interface between the outer and the inner parts 
of neutron star crust where unbound neutrons appear. At densities below ~ 10 12 
g.cm~ 3 , these unbound neutrons are not expected to be superfluid. 10 We will there- 
fore treat the nucleons in the Hartree-Fock approximation with the effective Skyrme 
nucleon-nucleon interaction. 11 Let us however point out that the above considera- 
tions are very general and apply to any approximation of the many-body problem. 
Indeed the Bloch states form a complete set of single particle states into which the 
many-body wave function can be expanded (for instance, the application of band 
theory including pairing correlations has been discussed in Ref. 12). 

The single particle wave functions are obtained by solving the self-consistent 
equations (q = n,p for neutrons and protons respectively) inside the Wigner-Seitz 
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where the single particle Hamiltonian is defined by 

h o ] =-V--^--V + U q (r)-iW q (r)-Vxa. (5) 



The effective masses m®(r), mean fields U q (r) and spin-orbit terms Wq(r) depend 
on the occupied single particle wave functions. 

It follows from the Floquet-Bloch theorem (fl]), that the wavefunctions can be 
expressed as 

^(r) = e* r ^(r), (6) 

where u^ q )(r) has the full periodicity of the lattice, u l(r + T) = u 1{t). Conse- 
quently, Eqs Q can be equivalently written as 

(^+^ ) )«S2W=e2«2(r) (7) 
where the fc-dependent Hamiltonian is defined by 

and the velocity operator v q is defined by the commutator 

v q ^±[r,h^}. (9) 



The boundary conditions are completely determined by the assumed crystal sym- 
metry. According to the Floquet-Bloch theorem fl} , the wavefunctions tp 1 (r) and 
Vak ( r + between two opposite faces of the Wigner-Seitz cell (T being the cor- 
responding lattice vector) are shifted by a factor e lkT while u^(r) = U^l(r +T). 
Each cell is electrically neutral which means that it contains as many protons as 
electrons. Equations are solved for a given baryon density, assuming that the matter 
is in /3-equilibrium. 

3. Validity of the Wigner-Seitz approximation 

The three dimensional partial differential Eqs ([JJ or (0 have to be solved for each 
wave vector k inside one irreducible domain of the first Brillouin ZoneQ. Such cal- 
culations are computationally very expensive (numerical methods that are applica- 
ble to neutron star crust have been discussed by Chamel 13 ' 14 ). Since the work of 
Negele&Vautherin, 5 the usual approach has been to apply the Wigner-Seitz approx- 
imation. 6 The complicated cell is replaced by a sphere of equal volume as shown in 



a The Hartree-Fock equations have to be solved together with Poisson's equation for determining 
the Coulomb part of the proton mean field potential. 
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Fig. [T] In this approximation, the dependence on the crystal structure is therefore 
completely lost. Besides only the states with k = are considered or in other words 
the Hamiltonian is neglected. Eqs (UJ) and |(7J) are then similar and reduce to a 
set or ordinary differential equations. The price to be paid for such a simplification 
is that this approximation does not allow the study of transport properties (which 
depend on the fc-dependence of the states). As pointed out by Bonche&Vautherin, 15 
two types of Dirichlct- Neumann mixed boundary conditions are physically plausible 
yielding a more or less constant neutron density outside the cluster: either the wave 
function or its radial derivative vanishes at the cell edge, depending on its parity. 

In recent calculations, 16-18 the Wigner-Seitz cell has been replaced by a cube 
with periodic boundary conditions. While such calculations allow for possible de- 
formations of the nuclear clusters, the lattice periodicity is still not properly taken 
into account since the fe-dependence of the states is ignored. Besides the Wigner- 
Seitz cell is only cubic for a simple cubic lattice and it is very unlikely that the 
equilibrium structure of the crust is of this type (the structure of the crust is ex- 
pected to be a body centered cubic lattice 1 ). It is therefore not clear whether these 
calculations, which require much more computational time than those carried out 
in the spherical approximation, are more realistic. This point should be clarified in 
future work by a detailed comparison with the full band theory. 

It should be remarked that the states obtained in the Wigner-Seitz approxi- 
mation do not coincide with those obtained in the full band theory at k = 0. 8 
The reason is that in the Wigner-Seitz approximation the shape of the exact cell 
is approximated by a sphere. By neglecting the fe-dependence of the states, the 
Wigner-Seitz approximation overestimates the neutron shell effects and leads to 
unphysical fluctuations of the neutron density, as discussed in details by Chamel et 
al. s Negele&Vautherin 5 proposed to average the neutron density in the vicinity of 
the cell edge in order to remove these fluctuations. However it is not a priori guar- 
anteed that such ad hoc procedure did remove all the spurious contributions to the 
total energy. This issue is particularly important since the total energy difference 
between different configurations may be very small. It has thus been found that the 
dependence of the equilibrium crust structure on the choice of boundary conditions, 
increases with density. 7 ' 15 In particular, the Wigner-Seitz approximation becomes 
very unreliable in the bottom layers of the crust where the clusters nearly touch. 

While the Wigner-Seitz approximation is reasonable at not too high densities for 
determining the equilibrium crust structure, the full band theory is indispensable 
for studying transport properties. 

4. Band theory and transport properties 

At low temperature the transport properties of the neutron liquid in neutron star 
crust are determined by the shape of the Fermi surface, defined by the locus of points 
in fc-space such that e a k = £f (we drop the superscript n since we only consider 
neutrons in this section). In general, the Fermi surface is composed of several sheets 
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Fig. 2. Neutron Fermi surface in the inner crust of neutron star at average mass density p = 
7 X 10 11 g.cm - 3 . The crust is composed of a body centered cubic lattice of zirconium-like clusters 
(Z = 40) with N = 160 neutrons per lattice site (70 neutrons are unbound). The single particle 
states have been calculated by solving the Hartree-Fock equations with the Skyrme force SLy4 
and by applying Bloch boundary conditions. 8 The Fermi surface is shown inside the first Brillouin 
zone. The different colors correspond to different bands. 

corresponding to the different bands that intersect the Fermi levei. An example of 
neutron Fermi surface is shown on Fig. O 

The motion of the unbound neutrons is affected by the presence of the nuclear 
clusters. In particular, they can be Bragg-reflected by the crystal lattice. As a result 
of the momentum transfer with the lattice, a neutron in a state characterized by a 
wave vector A; and a band a, move in the crystal as if its mass was replaced by an 
effective mass given by (considering cubic crystals) 



where denotes the Laplacian operator in fc-space. This concept of effective mass 
has been very useful in interpreting neutron diffraction experiments. 19 In the context 
of neutron star crust, the number of unbound neutrons can be very large and it is 
therefore more appropriate to introduce a "macroscopic" effective mass m* averaged 
over all occupied single particle states 



where rif is the number density of free neutrons and the integral is performed over 
all occupied states "inside" the Fermi surfac^). Let us remark that the mobility 
coefficient JC can be equivalently written as 



nin\ak = 3ft 2 (Afe£ afc ) 



-l 



(10) 




(11) 




(12) 



b The occupied states are those for which the energy e a k is lower than the Fermi energy ep. 
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showing that the effective mass depends on the shape of the Fermi surface. This 
effective mass to* has been found to be very large compared to the bare mass. 13 ' 14 
For instance, the effective mass corresponding to the Fermi surface shown in Fig. [3 
is about m* ~ 4.3 m n . On the contrary, this effective mass is slightly lower than the 
bare mass in the liquid core owing to the absence of clusters 20 (the effectice mass 
m* still differ from m n due to neutron- neutron and neutron- proton interactions). 
With the definition pi[) . it can be shown that in the crust frame the momentum 
Pi of the neutron liquid is simply given by pf — m^Vf, where Vf is the neutron 
velocity. 21,22 This result is quite remarkable. It implies that in an arbitrary frame, 
the momentum and the velocity of the neutron liquid are not aligned! Indeed, if v c 
is the velocity of the crust, applying the Galilean transformation leads to 

Pi = to* v f + (m n - m* )u c . (13) 

These so-called entrainment effects are very important for the dynamics of the 
neutron liquid in the crust. For instance, it has been suggested that for sufficiently 
large effective mass m* 3> m n , a Kelvin- Hclmholtz instability could occur and 
might explain the origin of pulsar glitches. 23 

5. Conclusion 

Since the work of Negele&Vautherin, 5 the Wigner-Seitz approximation has been 
widely applied to study neutron star crust. Nevertheless, the necessity to go beyond 
has become clear in the last few years. It has thus been shown that the results of 
the Wigner-Seitz approximation become less and less reliable with increasing den- 
sity. 7 Besides this approach, which decomposes the crust into a set of independent 
spherical cells, does not allow the study of transport properties. A more realis- 
tic description of the crust requires the application of the band theory of solids. 8 
By taking consistently into account both the nuclear clusters which form a solid 
lattice and the neutron liquid, this theory provides a unified scheme for studying 
the properties of neutron star crust. We have shown in this framework that the 
unbound neutrons move in the crust as if they had an effective mass much larger 
than the bare mass. The analogy between neutron star crust and periodic systems 
in condensed matter suggests that some layers of the crust might have properties 
similar to those of band gap materials like semiconductors or photonic crystals for 
instance. The existence of such "neutronic" crystals in neutron star crusts opens a 
new interdisciplinary field of research, at the crossroad between nuclear and solid 
state physics. 
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